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Introduction 

This  review  Is  confined  to  the  theory  of  laminar  shocks  in 
where  there  is  no  reference  to  instabilities  or  turbulence. 
This  restriction,  plus  the  limitation  to  perpendicular 
shocks,  limits  the  field  to  a  small  part  of  plasma  shocks. 
Nevertheless,  this  restricted  subject  is  sufficiently 
ramified  to  offer  a  very  wide  choice  of  physical  mechanisms 
and  a  theory  that  is  far  from  complete. 

The  original  colllslonless  shock  theories—'  were  laminar 
for  the  good  reason  that  they  antedated  all  experiments; 
it  is   a  very  appropriate  use  of  Occam's  razor  to  dismiss 
turbulence  until  it  is  forced  upon  one.'   Although  instability 
and  turbulence  are  observed  in  many  experiments,  they  are 
not  necessary  for  irreversible  colllslonless  heating,  and 
it  is  very  likely  that  nonturbulent  mechanisms  will  ultimately 
play  a  role  in  the  subject. 

The  key  point  is  that  there  is  an  enormous  selection  of 
possible  Irreversible  mechanisms  within  each  of  the  broad 
categories — collislonal,  laminar  colllslonless,  and  turbulent. 
Even  the  classical  collislonal  approach  offers  so  large  a 
variety  of  mechanisms  that  the  theory  is  far  from  complete. 
An  exception  is  the  case  of  weak  shocks  where  the  theory  is 
relatively  complete  and  the  complex  interplay  of  collislonal 
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mechanisms  can  be  viewed  explicitly.   Turbulent  theories 
of  shock  structure  can  be  expected  to  be  as  diversified  as 
plasma  stability  theories.   It  Is  certain  that  many  shocks 
will  exhibit  a  combination  of  two  or  more  mechanisms. 

The  significance  of  shock  waves  for  heating  plasma  Is 
that  the  plasma  Itself  decides  which  Is  the  most  effective 
mechanism.   Except  for  subtleties  of  bifurcation  and  non- 
unique  profiles  (of  which  there  Is,  so  far,  no  evidence)  the 
dominant  mechanism  will  be  the  one  which  gives  the  widest 
profile.   Roughly  speaking,  a  theory  which  predicts  too  wide 
a  profile  Is  In  error,  and  one  which  predicts  too  narrow  a 
profile  Is  Incomplete. 

We  give  a  brief  survey  of  the  subject  In  terms  of  the 
relevant  dlmenslonless  parameters;  then  an  account  of  the 
highly  developed  weak  shock  colllslonal  structure,  and  of 
the  scattered  results  which  are  known  In  colllslonless  theory. 
With  regard  to  dlmenslonless  parameters,  the  main  result 
Is  that  there  are  very  many.   With  regard  to  the  colllslonal 
theory,  a  significant  point  Is  that  experiments  are  only 
marginally  colllslonless,  by  a  factor  of  about  ten;  the 
theoretically  estimated  laminar  thickness  may  be  In  error 
by  a  comparable  amount.   The  contribution  of  laminar 
colllslonless  theory  Is,  at  present,  confined  to  substructures 
such  as  the  "foot"  ahead  of  a  sufficiently  strong  shock. 
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Parameters 

Counting  Is  the  most  basic  of  mathematical  operations; 
counting  parameters  In  a  physical  problem  Is  the  first,  and 
sometimes  even  the  most  important  mathematical  theory. 
Karl  Schlndler  described  the  problem  of  the  solar  wind 

flowing  past  a  spherical  earth  as  one  with  nine  dlmenslonless 

2/ 
parameters.—'   A  steady.  Infinite,  plane  shock  wave  Is  a 

simpler  physical  problem  with  no  physical  boundaries  and  no 

externally  Imposed  length  or  time  scales — It  therefore 

yields  "only"  six  dlmenslonless  parameters  (In  the  simplest 

version  of  a  fully  Ionized,  two  component  plasma,  with  unit 

Ionic  charge). 

From  the  molecular  constants, e  and  m,  the  electro- 

o 
magnetic  parameters,  f^   \i     =   l/c  ,  the  thermodynamic  state, n 

and  kT,  and  the  field  strength  H  (six  dimensional  parameters 

In  all),  we  can  construct  one  natural  mass,  length,  and  time, 

and  three  further  dlmenslonless  combinations.    We  can  take 

o 

these   to  be    P  =^  Sirp/H   ,    wt,    and  "V    =  <*)  /a>  .      To  account   for 

P  c 

o 

Ions  and  electrons,  we  add  a  =  m_/m  .   These  four  dlmen- 
slonless parameters  describe  a  state  of  thermodynamic 
equilibrium.   For  the  shock  we  add  a  strength  parameter  and 
an  obliqueness  parameter  to  make  six. 

In  a  perpendicular  shock  we  count  five  Instead  of  six. 
In  a  collislonless  shock  we  omit  oyr ;    in  most  colllslonal 
theories  -7  does  not  appear — therefore  four  parameters  is 
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common.  But  even  the  "small"  number,  four,  represents  a 
very  high  dimensional  parameter  space,  and  we  can  expect 
qualitatively  different  modes  of  behavior  to  occur  as  we 
vary  these  parameters. 
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Weak  Collislonal  Shocks 

The  one  case  in  which  there  is  a  relatively  complete 
theory  of  collislonal  shocks  is  that  of  a  weak  shock.   In  this 
limiting  case  an  explicit   solution  has  been  obtained  for 
the  system  of  two  Boltzmann  (or  Fokker-Planck)  equations 
and  Maxwell's  equations.—'   All  parameters  other  than 
shock  strength  are  arbitrary. 

For  the  Navier-Stokes  equations  in  an  ordinary  gas 
the  weak  shock  profile  is  given  as—'—' 

0-^6  tanh  (  ex/L)  (1) 

where  ^  is  density,  temperature,  or  velocity,  (suitably 
normalized)  and  e  is  the  shock  strength.   The  scale  length, 
L,  is  a  sum  of  two  mean  free  paths,  one  based  on  viscosity 
and  the  other  on  heat  conductivity. 

In  the  plasma  problem,  the  weak  shock  profile  is  given 

by  the  same  formula  (1),  except  that  the  scale  length  L  is 

"5/ 
a  very  complicated  algebraic  expression.—"^   This  expression 

is  greatly  simplified  by  taking  a  singly  charged  ionic 

species  and  letting  a  —*■  0,—' 
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where 


1 


^-   /208  ^  192  ,  3ajo^   ,  13  ?  ,2^2\/^ 
5408  ^  128  2^2  _^  ^2  4^4 


(3) 


and 


2  2 
3   k.  1 


m 


4l2¥  ne  log  Z 


(4) 


Is  a  standard  mean  free  path. 

Still  further  simplification  is  obtained  in  the  limit 
P  — f-   0,  ooT — >  00;  in  terms  of 

X  =  aooT 

y    =    P<OT 

we  have 


(5) 


T/T  2     n,f  6    ^^i2  4+25x^/181  (gs 

VLm==9  ^^  I  xy-^T8-, ^23x2/18^ 


-6- 


The  reason  for  the  complex  formula  (2)- (3)  lies  in 
the  complexity  of  the  physical  mechanism.   There  is  a  2  x  2 
matrix  of  viscosity  coefficients  coupling  ions  and  electrons; 
there  is  a  3  x  3  matrix  of  heat  conductivity,  diffusion 
(i.e.  resistivity),  and  thermal  diffusion  coefficients; 
in  addition  there  are  distinct  Hall  coefficients  for  ion 
and  electron  stresses  and  for  ion  and  electron  heat  flows 
as  well  as  the  usual  one  for  electric  current.   There  are 
more  than  a  dozen  distinct  regions  in  theO,t»>T)  parameter 
plane  where  different  dissipatlve  mechanisms  dominate.—'^ 

Figure  1  displays  contours,  L/L  =  const.,  for  P 
ranging  from  .001  to  10  and  u>r   from  1  to  1000  (a^  =  I/3672 
for  Deuterium).   The  scale  length,  L,  ranges  from  one- 
hundredth  to  one  thousand  times  the  mean  free  path  in  this 
region.   A  "standard"  calculation  of  the  shock  profile, 
using  simple  resistivity  alone,  gives  the  scale  length 


L3=  C^/ilTTCU  =  ^^  1^  (7) 


Figure  2  displays  contours  of  L/L  for  the  same  range  of  P 

and  COT  as  in  Figure  1.   The  value  of  L/L^  (ranging  from  one  to  ten 

thousand)  can  be  interpreted  as  the  value  of  an  effective 

collision  frequency,  v  f>f/v  ,  required  to  correct  the 

"standard"  collisional  thickness  L  .   The  standard  formula 

(7)  is  frequently  used  to  reduce  experimental  data,  giving 

a  value  of  v  „„. 
ef  f 

The  monotone  profile  (1)  results  if  e  is  sufficiently 
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small.   But  this  limit  Is  not  uniform  In  parts  of  theOjo:)^)- 
plane,-^  specifically.  If,  x  =  acox  becomes  large  and  y  =  ^'c 
becomes  small  as  e  — >0.   Instead  of  one  parameter  L,  the 
profile  depends  on  L  and  on 

A  =  (2€/5)V2  Lp^  Lp  =  cA)p^  (8) 

2   2 
For  M-  =  A  /L  <  1/8  the  profile  is  monotone;  but  it  displays 

damped  oscillations  when  [x   >   l/8  or 

Presumably,  in  a  time  dependent, steepening  wave  the  profile 
is  determined  by  the  first  (i.e.  the  larger)  of  the  two  length 
scales  L  and  X.    If  dispersion,  given  by  A,  interferes  before 
dissipation,  measured  by  L,  takes  over,  the  profile  oscillates. 
The  dashed  lines  in  Figure  1  mark  the  regions  where  oscillatory 
profiles  make  their  appearance  for  various  values  of  e. 

Another  nonuniform  limit  occurs  in  the  same  part  of 
the  parameter  plane  (viz.  x-^oo  and  y-*0)  when  -y  =  oj  /co 
becomes  small  with  e.   Another  oscillatory  profile  occurs, 
but  this  has  not  yet  been  calculated  quantitatively.—' 

Collisional  shock  profiles  for  finite  shock  strength 
have  only  been  calculated  for  a  limited  range  of  other  para- 
meters and  with  greatly  simplified  dissipative  mechanisms. ■^''^ 
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Since  the  more  sophisticated  dissipative  mechanisms  may- 
give  rather  large  values  of  v  ^^/v    ,    the  decision  that  a 
given  experiment  is  colllsionless  should  never  be  made 
solely  on  the  basis  of  a  "standard"  resistive  thickness 
or  an  equivalently  crude  estimate  of  the  collisional 
dissipation. 
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Colllslonless  Shocks 

The  universal  mechanism  for  all  Irreversible  behavior 
Is  that  small  changes  In  Initial  state  become  large  at  some 
later  time.   An  ordinary  molecular  collision  Is  an  extremely 
effective  mechanism  for  losing  Information  about  the  state 
of  an  Individual  particle.   Turbulent  electric  fields  can 
accomplish  a  similar  purpose  In  "colllslonless"  shocks. 
But  neither  collisions  nor  turbulence  are  necessary  for 
Irreversibility.   The  temi  phase  mixing  covers  a  multitude 
of  such  effects.   In  a  linear  system,  colllslonless  damping 
Is  evidenced  by  the  presence  of  a  continuous  spectrum--e.g. 
as  In  Landau  damping.   The  case  of  flow  perpendicular  to 
the  magnetic  field  Is  more  subtle.   Since  the  linearized 
equations  exhibit  a  discrete  spectrum,  dissipation  requires 
a  nonlinear  analysis.   Also,  in  the  transverse  steady  shock 
problem  it  can  be  shown  that  two  particles  which  are 
initially  close  will  never  be  separated  by  more  than  a 
finite  distance.   Nevertheless,  two  particles  which  are 
initially  close  can  suffer  entirely  different  histories  if, 
e.g.,  one  is  reflected  by  a  potential  and  the  other  is  not; 
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(the  reflected  particle  must  eventually  penetrate  the 
barrier  after  a  second  or  further  return).   This  mechanism, 
enhanced  by  the  presence  of  large  electric  fields  as  a 
consequence  of  a  small  mass  ratio  expansion,  has  been  shown 
to  give  rise  to  a  definite  irreversible  profile,  although 
not  a  complete  shock — specifically  through  a  transition 
from  a  constant  state  to  a  periodic  state  of  different 
entropy. — '       The  "foot"  which  has  been  observed  to  develop 
in  front  of  many  medium  strength  shocks  was  also  theoretically 
predicted  by  this  analysis. — ■ 

One  theoretical  avenue  for  an  irreversible  transition 
is  through  a  transfer  from  one  branch  to  another  of  a 

multivalued  direction  field. — '    There  is  precise  agreement 

12/ 
in  such  a  case  between  theory  and  numerical  experiments. — ' 

"Critical"  Mach  numbers  have  been  observed  experimentally 
across  which  a  qualitative  change  in  structure  is  observed. 
There  are  also  several  such  theoretical  transitions,  e.g. 
by  a  change  in  the  nature  of  a  singularity,  but  there  is  no 
definite  correlation  of  theoretical  and  experimental  transi- 
tions as  yet. 

Although  no  colllsionless  theoretical  model  has  yet 
given  a  complete  transition  from  one  constant  state  to 
another,  it  is  very  likely  that  they  will  play  a  role  in 
the  multlstructured  shocks  that  are  experimentally  observed. 
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A  Model  Equation 


An  equation  which  displays  many  of  the  qualitative 

13/ 
features  found  in  most  shock  waves  is— =^' 


U,  +mj=U+U  (10) 

t        X      XX      XXX 

All  of  the  coefficients  can  be  set  equal  to  unity  by 
scaling  x  and  t  and  U;  (but  this  merely  transfers  the 
criteria  distinguishing  different  solutions  to  the  initial 
data  ). 

The  two  terms  U,  +  UU  =0  exhibit  steepening,  U,  =  U 
damping.   The  combination  UU  =  U   gives  a  hyperbolic  tangent 
shock  as  in  (1),  and  the  combination  of  three  terms 
U^  +  UU  =  U   is  the  Burgers  equation  which  demonstrates 

t       X     XX  t!       -^ 


the  evolution  of  a  monotone  shock. 


ii/ 


The  pair  U,  =  U    demonstrates  simple  dispersion, 

UU  =  U    gives  periodic  and  solitary  wave  solutions,  and 
X     XXX  &      J" 

U,  +  UU  =  U    is  the  Korteweg  de  Vries  equation  which 

t       X     XXX  ^ 

describes  nonlinear  but  reversible  interaction  between 

15/ 
waves. -^^^"^ 

The  complete  steady  state,  UU^  =  U   +  U    is  the 

equation  mentioned  previously  which  describes  a  family  of 

steady  shock  profiles  ranging  from  monotone  to  damped 

oscillatory.   One  question  that  could  be  answered  by  a  study 

of  (10)  is  the  interplay  between  dissipation  and  dispersion 

in  a  steepening  wave,  as  previously  mentioned. 
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The  Initial  conditions  are  essential.   For  example, 
an  Initial  state  which  takes  the  same  constant  value  at 
X  :=  +  00  would  eventually  decay — '  ;    an  Initial  state  with 
different  constant  values  at  x  =  +  oo  should  develop  Into 
a  shock.   It  would  be  Interesting  to  study  the  evolution  of 
solutions  of  the  Korteweg-de  Vrles  equation  for  an  Initial 
state  joining  two  different  constant  states. 
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Fig.  1   Shock  Thickness  as  a  function  of  P  and  u)t  for  Deuterium 
(a^  =  1/3672) 
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